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Abstract. Power-suppressed corrections arising from end-point integration regions to the space- 
like vertex function of the massive ?7'-meson virtual gluon transition rj' — g*g* are computed. 
Calculations are performed within the standard hard-scattering approach (HSA) and the running 
coupling method supplemented by the infrared renormalon calculus. Contributions to the vertex 
function from the quark and gluon contents of the 7;'-meson are taken into account and the 
Borel resummed expressions for F^/g.g* (Q^, lj, 77), as well as for F^rggt{Q^,uJ = ±1,77) and 
Fjii g* {Q^ , ui = 0,ri) are obtained. It is demonstrated that the power-suppressed corrections 
~ (/IVQ^)", in the explored range of the total gluon virtuality 1 < < 25 GeV^ , considerably 
enhance the vertex function relative to the results found in the framework of the standard HSA 
with a fixed coupling. Modifications generated by the Tj'-meson mass effects are discussed. 

PACS. 12.38. Bx, 14.40.Aq, ll.lO.Hi 
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1 Introduction 

Recently interest in theoretical investigations of the 
gluonic structure of the 77 and rj' -mesons has risen 
due to the high precision CLEO results on the elec- 
tromagnetic ?77, T]'^ transition form factors (FFs) [1], 
as well as because of the observed large branching ra- 
tios for the exclusive B ^ K + rj' and semi-inclusive 
B^r]' +X, decays [2l|3]. 

The data on FF of the rj'j transition were mainly 
used for extracting constraints on the quark compo- 
nent of the ?7'-meson distribution amplitude (DA) [H 
[5l[6] . In these investigations various theoretical schemes 
and methods were employed. An important conclu- 
sion drawn from these studies is that the quark com- 
ponent of the Ty'-meson DA should be close to its 
asymptotic form and that the admixture of the first 
non-asymptotic term should be within the range of 
5|(1 GeV^) ~ 0.05 - 0.15, B| being the first Gegen- 
bauer coefficient. 
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An effect of the gluon component of the ry'-meson 
DA on the t]'^ transition was analyzed in [Till] , where 
relevant constraints on the input parameters i3| and 
B2 were extracted: it was shown that their allowed 
values are strongly correlated. Useful bounds on the 
Gegenbauer coefficients i?| and i?f were obtained 
also from investigation of the semi-inclusive decay 
r(15)^ry' + X fg]. 

The two-gluon valence Fock component of the rj'- 
meson can directly contribute to the rj'^ transition 
FF only at the next-to-leading order due to quark 
box diagrams and also affect the leading order re- 
sult through evolution of the quark component of 
the ?7'-meson DA. Hence, an effect of the r;'-meson 
gluon component on the r]'^ transition is mild. Con- 
trary, the contribution of the gluon content of the 
77'-meson to the two-body non-leptonic exclusive and 
semi-inclusive decay ratios of the _B-meson may be 
sizeable. Indeed, in order to explain the observed 
large branching ratio Bi{B rj' + Xg), in [10] a 
mechanism that employs the two-gluon content of the 
?7'-meson was suggested. In accordance with this ap- 
proach the dominant fraction of the B ^ ij' + 
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decay rate appears as the result of the transition 
9* ~* dv' of ^ virtual gluon from the standard model 
penguin diagram h — ^ sg* . In .lOj the g*gri' ver- 
tex function (VF) was approximated by the constant 
i/(g2,0,m2,) ~ 77(0, 0,m2,) ~ 1.8 GeV"\ the latter 
being extracted from the analysis of the J /if) -q'^ 
decay. Further investigations demonstrated that the 
effects of the QCD running couphng as(q^) p4l, as 
well as the momentum dependence of the form factor 
H(q^, 0,171^,), properly taken into account, consider- 
ably reduce the contribution to Br{B rj' + Xg) of 
the mechanism under consideration [12j . To eliminate 
the discrepancy between theoretical predictions and 
the experimental data in [13j a gluon fusion mecha- 
nism was proposed. In accordance with the latter, the 
77' -meson is produced by the fusion of a gluon from 
the QCD penguin diagram b — > sg* with another one 
emitted by the light quark inside the S-meson. In this 
mechanism, the vertex function -fr^'g'g* (<Zi^, 9I, "i^^,) 
appears owing to the g*g* —^ 77' transition. Similar 
ideas form a basis for the computation of the branch- 
ing ratios of various two-body non-leptonic exclusive 
decay modes and transition FFs of the B meson [141 
fT5lfTBlfT7lfT8lfT9l[20]. 

Hence, the Ty'-meson virtual (on-shell) gluon tran- 
sition VF, F^ig»g<.{q1,q2,m^,), is the central ingredi- 
ent of the relevant analysis performed within pertur- 
bative QCD (pQCD) and it deserves further investi- 
gations. This VF was computed in numerous works 
f8l[2Tl[22l[23]. The space-hke massless ?fg*g* vertex 
within the standard hard-scattering approach 
(HSA) was considered in [8], where correct analysis 
of the normalization of the gluon component of the 
jy'-meson DA and that of the gluon projector onto 
a pseudoscalar meson state was performed. Power- 
suppressed corrections, arising from the end-point in- 
tegration regions a; ^ 0, 1, to the massless space- like 
ri'g*g* VF were found in 22J. In this work the stan- 
dard HSA and the running couphng (RC) method 
together with the infrared (IR) renormalon calculus 
were applied. The 77'-meson mass corrections to the 
rj' g* g* space- and time-like vertex in the standard 
HSA were calculated in [23]. 

In the present work, we extend the results ob- 
tained in [22] by taking into account the jy'-meson 
mass effects, which may be considerable. The RC 
method [24ll25] enables us to estimate power correc- 
tions coming from the end-point a; — > 0, 1 regions 
in the integrals determining the amplitude of the 
1]' g* g* transition. Indeed, in the framework of the 
standard HSA [26l[27ll28] . in order to calculate the 
amplitude of the process, one has to perform inte- 



grations over the longitudinal momentum fractions 
of the constituents of the meson. If one chooses the 
renormalization scale in the hard-scattering am- 
plitude Th of the corresponding partonic subprocess 
in such a way as to minimize higher-order corrections 
and allows the QCD coupling constant Q!s(/i|j) to run, 
then one encounters divergences arising from the end- 
point a; — > 0, 1 regions. The reason is that the scale 
, as a rule, is equal to the momentum squared 
of the hard virtual partons carrying the strong in- 
teractions in the subprocess' Feynman diagrams and 
depends, in general, on x (or x = 1 — x). Within 
the RC method this problem is resolved by applying 
the renormalization group equation and the IR renor- 
malon calculus (for a review see, [S^I^ V It turns out 
that such treatment allows us to evaluate power cor- 
rections to the physical quantity under consideration 

[iiziiiMim^. 

This paper is structured as follows. In Sect. H] we 
present the necessary information on the quark-gluon 
structure of the Ty'-meson, its DAs and the hard- 
scattering amplitudes of the relevant subprocesses. 
Sect. [3] is devoted to calculation of the quark compo- 
nent of the T]' g*g* vertex function. The contribution 
to the VF of the gluon content of the Ty'-meson is 
computed in Sect. UJ Section [5] contains our numeri- 
cal results and their analysis. In Sect. |S] we make our 
concluding remarks. 



2 Quark and gluon content of the 

r^'-meson and the ri'g*g* vertex 

The Fock state decomposition of the pseudoscalar 
P ~ 7], ?7'-mesons can be written in the following 
form 

\P) = \Pa) + in) + \Pc) + \Pa) , 

where \Pa) and \Pb) denote the P-meson light quarks, 
and |Pc) and \Pg) its charm and gluon components, 
respectively. 

The light-quark content of the P-meson can be 
described either in the SUf{3) octet-singlet or in the 
quark-flavor basis. In this paper we choose to work 
in the quark-flavor basis 



Vq) ^ ^\uu + dd) , \r]s)=^s\s) 



(1) 



Here i'i denote wave functions of the corresponding 
parton states. 

We neglect the charm component of the ry'-nieson, 
because in accordance with existing estimations (34] 
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[35ll36j . it is too small to affect considerably the B- 
meson exclusive decays. 

The pure light-quark sector of the rj — rj' system 
without charm and gluon admixtures can be treated 
as superpositions of the basic states 



i't'p \Vq) - sin0p \r]s) 



\r]') ^ sin (l)p Ifjq) + cos (l)p\r]s) . (2) 

One of the advantages of the quark-flavor basis is 
that in this basis the decay constants /p"'' follow 
with great accuracy the pattern of the state mixing 

m 

fri = fq COS (t>p, =-fs Sm (t>p, 



fn' = f<} Sin'/'p, f^y = fs C0S(/>p, 



(3) 



where the decay constants fq and /s, and the mixing 
angle 4>p have the values 

/, = (1.07 ± 0.02)/,, /, = (1.34 ± 0.06)/,, 

= 39.3° ± 1.0°, (4) 

with /, — 0.131 GeV being the pion weak decay con- 
stant. 

The singlet part of the 77'-meson DA, which is 
only relevant to our present investigations, depends 
on both the quark (p''{x, fi^) and gluon 0^(x, /i^) com- 
ponents of the 77'-meson DA. These functions satisfy 
the symmetry and antisymmetry conditions under 
the exchange x ^ x, 

0«(x,M') = 0n^,M'), = (5) 

and they are given by the expressions 



09(x,/i^) = 6Cxx< 1 + J2 



^ rn n 



Q^s(a^o) 
as(M|) 



Cl'\x-x) 



(6) 



and 



n=2,4.. 



Qs(a*o) 
as(A*F) 



as(Atg) 



(7) 



where the constant C is defined as 

C = V2fq sin (j)p + fs cos 4>p. 

In dg and 0, Cf/^(z) and C^{.^i(z) are Gegen- 
bauer polynomials, fi^ and /Xg are the factorization 
and normalization scales, respectively. The values of 
the input parameters B"^ and B^ have to be fixed 
at the normalization scale (Iq — 1 GeV^: they deter- 
mine the shape of the DAs. In the above expressions, 
Q;s(/i^) is the QCD coupling constant in the two-loop 
approximation given by 



4n 



1 



2/3i lnln(/iVyl2) 

1[ Hf^yA^) 



(8) 



with Pq and Pi being the one- and two-loop coeffi- 
cients of the QCD beta function 

2 19 

f3o = n^-nf, p, = 51--nf. (9) 

Here, A is the QCD scale parameter and Uf is number 
of active quark flavors. 

In this work, we shall use the Ty'-meson DA that 
contains only the first non-asymptotic terms. Stated 
differently, we suppose that in ([6]) and (H]) i3| ^ 
0, S| ^ and B« = Bff = for all n > 4. Taking 
into account the expressions for the required Gegen- 
bauer polynomials as well as the values of the rele- 
vant parameters, we can recast the ?7'-meson quark 
and gluon DAs into the following simple forms [55] 

(j)'^{x, n%) = eCxx [l + A{^p) - 5A{fi%)xx] , 



4,a(x,fil) = Cxx{x-x)B{fil). (10) 

For Uf — 4 the functions A{^'^) and S(/ip) are de- 
fined by 



&Bl 



19B| 



Qs(m|) 

as(Mo) 

«s(a^f) 
as(Mo) 



B'' 
17 



bB?, 



as(MF) 
as(Mo) 



107 
75 



(11) 
3 can 



The ry'-meson quark and gluon DAs for Uf 
be found in [22] . 

The ?7'-meson virtual gluon transition vertex func- 
tion Fr^> g. {Q^ , uj , 1]) is the sum of the quark and 
gluon components defined in terms of the invariant 
amplitudes for the process 
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in the following way: 



41.1/prj a* b* 



(12) 

Here, e°, and qi, q2 are the polarization vectors 
and four- momenta of the two gluons, respectively. Be- 
cause we study only the space-like VF, qf and g| obey 
the constraints Qi — —qf > and Q2 = — (/I ^ 0- 
The VF, Fr^igtgt {Q^, ij), depends on the total gluon 
virtuality Q^, the asymmetry parameter to, and 77' - 
meson scaled mass 77 




7P Y 




Fig. 1. Leading-order Feynman diagrams contributing to 



— (13) the hard-scattering subprocess g + g <?* + g* 



The parameter w varies in the region — 1 < < 1- 
The values w = ±1 corresponds to the Ty'-meson on- 
shell gluon transition and w = to the situation when 
the gluons have equal virtualities. 

In accordance with the factorization theorems of 
pQCD, at high momentum transfer the components 
of the VF F^^gig-^ {Q^, ^, v) can be calculated by means 
of the formulas 



pi 
v'g'g 



Jo 

+T^{x, Q^, u, r,, fi^)] 0«(x, fil)dx, (14) 



and 



-Tl{x, Q\ uj, r,, /i|,)] 0f (x, ^il)dx. (15) 
The sum 

and the difference 

-Ti{x,Q'^,uj,T],^ip) 

are the hard-scattering amplitudes of the subprocesses 
q + q g* +g* and g + g g* +g*, respectively. The 
Feynman diagrams contributing at the leading order 
to these subprocesses are depicted in Figs. [Hand O 

At the leading order of pQCD, the hard-scattering 
amplitudes do not depend on the factorization scale 

but depend implicitly on the renormalization scale 





Fig. 2. Feynman diagrams contributing at leading order 
to the subprocess g + g ^ g* + g* ■ 



fiji through asifJqi). scales and fj,j^ are inde- 

pendent of each other and can be chosen separately, 
we adopt in this work the standard choice for the 
factorization scale fij^ = Q^, and we omit in what 
follows the dependence of the hard-scattering ampli- 
tudes on Thus, we have 



Ti{x,Q'^,u;,r],fi%) = - 



3g2 ujX 



uj{l + A) — i]{x — x) 
x{l + Lo) + x(l — Lu) + 2xxri ' 



(16) 



and 



Tf(x,Q^,tj,77,/i^) 



Q^uf 



x{l +Lo) + x{l -uj) + 2(1 + xx)r] 
LljX [x{1 + lu) + x{l — Lo) + 2xxrj\ ' 



(17) 



where A = (l + 2r//a;^ +77^/0;^)^/^ [23]. The remaining 
two functions can be obtained from (|16p and (flTl) by 
means of the replacement x f^x. 
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In the standard HSA one sets the renormalization 
scale to be — Q^, fixing as(/x|j) with respect to 
x, and simphfying the calculation of the VF consid- 
erably. In this approach the functions (fT6|) and (fT7| 
possess the symmetry features 

T^2)^x,Q^,-uj,ri) = -T^^^^{x,Q^,uj,r]). (18) 

Using these features, as well as ^ and the symmetry 
properties of T^2)^ under x x exchange, one can 
prove that 

F:;fi^4Q',^,v) - Fff,^,iQ',-u,^). (19) 

The last equality is a manifestation of the Bose sym- 
metry of the process under discussion under vector 
particles-gluons exchange. 

Within the RC method the renormalization scale 
is chosen equal, as a rule, to the momentum 
squared \q'^\ of the virtual partons in the correspond- 
ing Feynman diagrams. For the massless yy'-meson on- 
shell gluon transition the scale is exactly equal to 

= Q^x. (20) 

Then upon x ^ x exchange, the argument of ofg (/!/?) 

in the functions T-f becomes equal to Jlj^ — Q'^x. In 
the general case, the absolute value of the square of 
the four-momenta of the virtual partons depends 
on the total gluon virtuality Q^, the asymmetry pa- 
rameter CO and the scaled mass term 77. However, to 
avoid problems related to the appearance of the pa- 
rameters Q^, cu and 77 in the argument of as, we shall 
use the renormalization scale (|20l) . Such a choice is 
justified from the physics point of view as well, be- 
cause namely the part ~ Q^x of the renormalization 
scale leads to the power corrections ~ (yl^/Q^)", n ~ 
1, 2, ... to the VF Fjji g* {Q^ , uj , 1]) , which we are go- 
ing to compute. 

In the present paper we adopt the symmetrized 
RC method, where as{Q'^x) and as{Q^x) are replaced 

by 

as{Q^x)+as{Q^x) 
2 

After this modification all symmetry properties of 
the hard-scattering amplitudes and the vertex func- 
tion remain valid within the RC method as well. The 
symmetrized RC method was successfully employed 
in the investigation of various exclusive processes [22l 

mm- 



3 Quark component of the vertex 
function F^^,^.^,(Q^ ^, 77) 

The expression for the quark component of the ver- 
tex function F^,g,g,{Q'^,uj,r]) can be computed only 
after resolving the problems of the soft end-point re- 
gions x — *■ 0,1. In fact, having inserted the explicit 
expressions of the hard-scattering amplitude and the 
quark component of the Ty'-meson DA into Eq. (I14p . 
one encounters divergences, arising from the singular- 
ities of the coupling constants as{Q'^x) and as{Q^x) 
in the limits a; ^ 0, 1. The RC method provides the 
required prescription to cure these divergences. 

To this end, we express the running coupling 
as{Q^x) in terms of as(Q^) and, as a result, obtain 
integrals that can be regularized and calculated us- 
ing the approach described in pill^ . Then the quark 
component of the VF is written as a perturbative se- 
ries in as{Q^) with factorially growing coefficients. 
The resummation of such a series is performed by 
means of a Borel transformation. Namely, one has to 
determine the Borel transform of the corresponding 
series and in order to get the resummed expression 
for the vertex function one has to invert the Borel 
transform. The Borel transform of the series with fac- 
torially growing coefficients contains infrared renor- 
malon poles located at the positive axis of the Borel 
plane; therefore the inverse Borel transformation suf- 
fers from the pole divergences. In other words, the 
Borel technique transforms the end-point divergences 
into the IR renormalon pole divergences of the inverse 
transformation. Then, the resummed expression can 
be extracted by computing the relevant integrals in 
the sense of the Cauchy principal value (29u30j . 

A useful way to bypass these intermediate oper- 
ations and directly obtain the Borel resummed ex- 
pressions is to introduce the following formula for 
asiQ^x) i: 

a,(g2x) = / duexp{-ut)R{u,t)x~'', (21) 
Po Jo 

where the function R(u, t) is defined as 

2/3i 

R{u, t) — 1 ^""(1 — 7 — hit — Inu). 

Po 

In the above, 7 ~ 0.577216 is the Euler constant and 
t = ln(gVyl2). 

Calculations of the quark component of the VF 
lead to the following result: 

F'^,g,g, (Q^ OJ, v) = F^{Q\ oj, r,) + F^Q^u;, rj), 

(22) 
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where 



FiiQ ,UJ,1]) = K^{UJ,1^) 

W Po 

/"OO 

/ due-''*R{u, t)B{2 - u, 2) [2F1 (1, 2; 4 - w; n) 
Jo 



+,F, {1,2 -u-A-u; n)] + rj) 



/ due-''*R{u, t)B{3 - u, 3) [2F1 (1, 3; 6 - w; n) 
Jo 

+2Fi{l,3-u;6~u;n)]. (23) 

2ry 



Here 

Ki(u;,?7) 

and 



V 

Ll!X{luX — — 77) ' 



Ll!X{luX — uj + rj) 



UJ + T] — UjX ' 



(24) 



r2 = 



2r] 



rj — UJ + ujX 



(25) 

Equation ([22]) with the functions F;f and F2 is the 
Borel resummed expression for the quark component 
of the rj' g* g* vertex function and contains power- 
suppressed corrections coming from the soft end-point 
regions [22] , 

Since under the replacement oj —uj the rela- 
tions Ki ^ K2 and ri ^ r2 hold, we get 



V g 9 



The components {Q"^ , oj , rj) and F2 {Q'^ , uj , rj) of 
F^ig^g, {Q'^,uj, rj) can be obtained from each other uti- 
lizing the transformation 



2F1 (a, 6; c; z) = (1 - z) ''2F1 [a,c- b; c: 



z- 1 



(26) 

The argument ri of the hypergeometric functions in 
Fi{Q^,u}, rj) in the region ui S (—1, 0) obeys the con- 
straint ri < 1, whereas r2 < 1 in the domain lo £ 
[0,1). In order to reveal the IR renormalon struc- 
ture of the Borel resummed vertex function, as well 
as to perform numerical computations, we have to 
expand the hypergeometric functions over ri or r2- 



We choose to work in the region uj G [0, 1), and we 
therefore employ the expression 

F^%.g.(g^c.,77) = 2F|(Q^c.,,7). (27) 

Due to the symmetry of F^i^tg* {Q^, w, if) under lo ^ 
—Lo, the pole structure of the vertex function and 
numerical results in the region uj g (—1,0) are the 
same. 

For the ?7'-meson on-shell gluon transition, we get 



F^^,gg,{Q\u = ±l,r^) 



167r2C[l + A(Q2)] 1 



1 + ?; 



poo 

X / due~''*R{u,t) [B{l,2~u) 
Jo 

807r2CA(g2) 1 



+3(2,1 ~u}] 



Po l + r? 



/•CO 

X / dMe-"*i?(u,t)[B(3,2-u) + B(2,3-u)]. 
Jo 

(28) 

In the case of gluons with equal virtualities, = 0, 
the F^,g,g,{Q'',uj = 0,77) can be obtained from (gT]) 
upon the substitutions K2 ~* Kq and r2 ^ tq, where 

Koiv) 



r7v/TT2A;(l + v/l + 2/77)' 



ro 



1 + Vl + 2/r;' 



(29) 



As we have mentioned earlier, the massless rj'- 
meson virtual gluon transition vertex function was 
computed in |22j. The predictions of this work for 
F^,g,g, {Q'^,u}, rj) should lead to the results of [22] in 
the limit of 77 ^ 0. To regain these results, it is nec- 
essary to expand the relevant functions over rj <^ 1, 
and only after that take the limit 77 ^ 0. Then, in 
the general case [uj 7^ 0, ±1) 



K2{uj,r]^Q) 



2uj 



1 



''2 



1 



Using the last expressions, it is not difficult to check 
that ([271) coincides with (4.4) of [22]. By setting 77 = 
0, from (|28|) one can easily recover the expression for 
the massless 7;'-meson on-shell gluon transition VF 
derived in (4.5) of [22] . 

The important question to be clarified here is 
whether one can use the results obtained within the 
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RC method in the hmit ^ oo in order to re- 
gain the asymptotic form of the VF. Indeed, regard- 



less of the method used, in the hmit Q 



the 



VF must reach its asymptotic form, because power- 
suppressed corrections vanish in the asymptotic hmit. 
In the hmit — > cxd, the gluon component of the 
T^'-meson DA vanishes, (j)^{x,Q'^) 0, whereas the 
quark component (j)''{x,Q'^) evolves to its asymptotic 
form 

0'?(a;,Q2) ^ QCxx. 

Therefore, the results that we obtain here not only 
describe the asymptotic limit of the quark component 
of the VF, but the asymptotic limit of the full VF 
itself. 

From the whole analysis performed in [22 , it fol- 
lows that in the asymptotic limit the substitution 



An 



due ^^R{u,t)B{n,m- u) 



X2-F'i (1, rn — li; TO + n — u; r) 

as{Q^)B{n, m)2Fi (1, m; m + n; r) , 

has to be appUed. 

Having used this prescription, we obtain 



X2i^i(l,2;4;r,). 



(30) 



Equation (|30p with the quantities labeled i = 1, 2, in 
general, should be employed in the relevant regions 
of the asymmetry parameter, i.e. in the regions uj e 
(—1,0) and U! G (0, 1), respectively. But, because the 
hypergeometric function 2^'i(l, 2, 4; 2;) is expressible 
in terms of the elementary ones, one can use ([5D|) with 
both i = 1 and i = 2 in the whole region uj E (—1, 1), 
excluding the point oj = 0. At uj ^ I {uj = —1) ([30]) 
with i = 2 [i = 1) can be applied. Our formula for the 
asymptotic limit of the quark component of the VF 
numerically is identical to (60) of [22 (after setting 
there y^/,,' C, A2{Q^) ~ and evolving the 
argument of to Q^). 

The IR renormalon structure of the expressions 
and (pS]) does not differ from that of the cor- 
responding massless r/'-meson gluon transition FFs 
described in rather detailed form in [52] , to which we 
refer interested readers. 



4 The gluon component of the vertex 
function 

We compute the gluon component of the VF, 
Ffi'g'g' ^1 ^ employing the formula 

^„^'<,-,*(Q',^,^) = 2 / Tf{x,Q\uj,r^)4><>{x,Q^)dx, 



which leads to the following result: 



(Q^ ry) = F^^{Q' , c., r;) + F^{Q\ u, rj). 

(31) 

In (|3ip the a and b components are given by the ex- 
pressions 



2 _ 4tt^CB{Q^^ 



V 



Q'^fion f (cjA)2 (tj + 77 + ljA) 



X dwe""*i?(u,t){(l + cj)[S(4-u,2) X 

2Fi(l, 2; 6-u;r) + B{4, 2 - u)2Fi (1, 2 - u; 6 - u; r)] 

-(1-w) [B(4,2-u)2Fi(l,4;6-«;r) 

+B{A - u, 2)2^1 (1, 4 - u; 6 - m; r)] 

-2u!B{3, 3-u) [2F1 (1, 3 - u; 6 - w; r) 

-f2Fi(l,3;6-u;r)] 

+277 [5(2, 3 -li) (2Fi(l,2;5-u;r) 

-2i^i (1, 3 - u; 5 - u; r)) -f S(3, 2 - u) 

X (2F1 (1, 2 - m; 5 - u-r) - 2F1 (1,3; 5 - u;r)) 

+5(4- w, 3) (2F1 (1,3;7-M;r) 

-2F1 (1, 4 - u; 7 - u; r)) + 5(4, 3 - u) 

X (2F1 (1, 3 - m; 7 - u; r) - 2F1 (1, 4; 7 - u-r))]} , 

(32) 

where 

277 



cj + 77 -I- wA ' 



and 



Q'^PoUf {ujX)^{r]-uj + ujX) 



X dMe-"*i?(u,t){(l +w) [5(4- u,2) X 
Jo 

2F1 (1, 4 - w; 6 - u; r) + 5(4, 2 - M)2i^i (1, 4; 6 - u; r) 
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-(1 - c^) [B(4, 2 - u)2Fi (1, 2 - It; 6 - u; r) 
+B(4-u,2)2Fi(l,2;6-u;r)] 
-2ujB{3, 3 - li) [2F1 (1, 3 - u; 6 - u; r) 
+2F1 (1, 3; 6- It; r)] + 277 [B(2,3-m)x 
(2F1 (1, 3 - u; 5 - u- r) - 2F1 (1, 2; 5 - u; r)) 
+B{3,2-u) (2i^i(l,3;5-it;r) 
-2F1 (1, 2 - u; 5 - It; r)) + 5(4 - u, 3) 
X (2F1 (1, 4 - u; 7 - m; r) - 2F1 (1,3; 7 - w; r)) 
+5(4, 3 -u) (2i^i(l,4;7-u;r) 

-2Fi(l,3-ti;7-u;r))]}, (33) 

with r = r2. 

It is worth noting that F^{Q^,uj,ri) has been ob- 
tained from the original result after the transforma- 
tion ([^5]) . In the region uj £ [0, 1) both r and r < 1, 
where (|32p and p3)) can be used for expansion and 
numerical calculations. We have checked that, 

which becomes evident if we represent the gluon com- 
ponent of the VF in the form 

F^, {Q\uj, r?) = Tf (x, Q\uj, rj)cj,Hx, Q^)dx 
Jo 

- [ Tl{x,Q\oj,r^W{x,Q^)dx. 
Jo 

But expressions obtained using Tj* (x, Q^, to, 77) are as 
lengthy as the ones presented in (|32p and (|33p : there- 
fore, we refrain from writing them down here. 

In the case of the Ty'-meson on-shell gluon tran- 
sition, the VF can be found after changing the fac- 
tors and arguments of the hypergeometric functions 
in (1211) and i.e. 
in FS{Q^,uj, 7^) 

V V - . V 

(a;A)2(w -I- wA) ^ 2(1 + 77)3 ' ^ l + 

mFl^{Q^,uj,7j) 

V ^ 1 ^ -, 

(WA)2(77 -LU + LUX) ^ 2(1+ ^ ■ 

The gluon component of the VF is identically equal 
to zero for equal gluon virtualities (w = 0) . This is the 



important qualitative modification induced by the rj'- 
meson mass term kept in the hard-scattering ampli- 
tudes. Let us emphasize that F^,g,g, (Q^, = 0, 77) = 
within both the standard HSA and the RC method. 

In the limit 77 — > our results reproduce the ex- 
pression for the gluon component of the massless 77'- 
meson virtual gluon transition VF from (4.22) of [22] . 
In fact, acting as in the case of the quark component 
of VF, we can show that the factor in F^ {Q'^ , uj , rj) 
vanishes, 

H .0 

(cjA)2(tj + 7/ + cjA) 

and for Fjl {Q'^ , ui , rj) get: 

?7 1 2w 

(a;A)2(77 — a; + a;A) a;(l+a;)' 1 + lo 

Then it is easy to demonstrate that, the function 
F^{Q'^,ui, 77) in the limit 77 ^ leads to F^,^,^, (Q^ uj). 

The infrared renormalon structure of the terms in 
([32| and p3|) ~ (1 -I- uj), (1 — uj), 2uj are the same 
as in the case of the massless 77'-meson virtual gluon 
transition VF. The terms ~ 277 are new; nevertheless, 
their IR renormalon structures can be clarified using 
the procedures described in 



5 Numerical analysis 

In order to start numerical computations, we need 
to fix the values of some constants and parameters. 
In our calculations the 77'-meson mass is set equal 
to 771^/ = 0.958 GeV. The value of the QCD scale 
parameter for 7^/ = 4 is vl = 0.3 GeV. 

To proceed with the computation of the 77'-meson 
gluon vertex function and explore the role played by 
the 7;'-meson gluon content and its mass in this pro- 
cess, we have to define also the allowed values of the 
input parameters i?| and i3| at the normalization 
scale fJ-Q — 1 GeV^. In the present paper we use the 
77'-meson asymptotic DA or select values of the pa- 
rameters B2 and B2 that obey the constraints 

51 = 0.1, Sfe [-2,14]. (34) 

The quark component of the 77'-meson virtual gluon 
transition VF for different values of the asymmetry 
parameter is shown in Fig. [3^. The chosen values 
of the input parameters correspond to the 77'-meson 
asymptotic DA. Since, in the case of the asymptotic 
DA, the gluon component of the VF vanishes, in this 
figure we, actually, have curves for the full VF. In the 
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Fig. 3. The quark (a) and gluon (b) components of the 
scaled VF Q'^ F,^i g* g, {Q'^ ,uj,ri) as functions of Q^. The 
solid curves are obtained using the RC method, whereas 
the dashed lines are calculated within the standard HSA. 
In b the correspondence between the curves and input 
parameter B| is _B| = 8 for the curves 1 and Bf = 4 for 
the curves 2 



Fig. 4. The quark (a) and gluon (b) components of the 
VF versus Q^. All curves are obtained in the context of 
the RC method. The solid lines are calculated by taking 
into account the Ty'-meson mass effects: in computations 
of the dashed lines the ry'-meson mass term is neglected. In 
b the correspondence between the curves and parameter 
B| is the same as in Fig. [3] 



same figure predictions obtained within the standard 
HSA are also depicted. One sees, that in the domain 
1 < Q2 < 25 GeV^ the standard pQCD results get 
enhanced by approximately a factor of two due to 
power corrections. A similar conclusion is valid also 
for the gluon component of the VF for > 4 GeV^ 
( Fig. [St) as weU. 

Even from these first results it is evident that 
soft end-point corrections lead, approximately, to the 
same enhancement of the standard predictions, as in 
the case of the massless j^'-meson virtual gluon tran- 
sition [22]. Therefore, it is interesting to find modi- 
fications in the behavior of the VF induced by the 
ry'-meson mass effects. 



In Fig.[31 the VFs computed within the RC method 
by taking into account and neglecting the ry'-meson 
mass term, are demonstrated. The differences in be- 
havior of the quark component (Fig. 4a) are con- 
siderable. Indeed, mass effects suppress the absolute 
value of the quark component and, at the same time, 
change its shape in the pQCD applicable region of 
Q^. The gluon component, as a function of the total 
gluon virtuality Q^, is not affected dramatically by 
the ?7'-meson mass effects (Fig. 4b). 

But the Ty'-meson mass term changes drastically 
the behavior of the gluon component of the VF as a 
function of the asymmetry parameter. As a function 
of ui, the gluon component is plotted in Fig. [Sb- It 
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Fig. 5. The quark (a) and gluon (b) components of the 
VF at fixed -B| and versus lo. In a. all curves are ob- 
tained in the context of the RC method. The correspon- 
dence between them and the parameter Bf is: _B| = for 
the solid curve; _B| = 4 for the dashed curve; _B| = 8 for 
the short-dashed curve. In b the solid curves are obtained 
within the RC method. For computation of the broken 
lines the standard HSA is used. The correspondence be- 
tween the curves and the parameter Bf is i?| = 8 for the 
curves 1; Bf — 4 for the curves 2 and Bf = 2 for the lines 



turns out that F^,g,g,{Q^ ,uj,f]) vanishes a.t lu — 
in the framework of both the standard HSA (the 
dashed hues) and the RC method (the sohd hnes). 
For small uj the end-point effects are also mild. There- 
fore, it is legitimate to conclude that in the region 
|w| < 0.2, F^, g, g, {Q"^ , oj , rj) feels neither end-point 
nor mass effects. The dependence of the quark com- 
ponent F^i gt g, (Cp T ^ T Tf) on LO is plotted in panel (a) 
of the same figure: as a function of lj it demonstrates 
firm stability against variations of i?| . 



Fig. 6. The quark (a) and gluon (b) components of the 
VF as functions of at lj = ±1. All curves are computed 
using the RC method. The correspondence between plot- 
ted curves and parameter B| is B| = for the solid 
curves; B| = 2 for the dot- dot- dashed curves; B| = 4 
for the dashed lines; Bf = 6 for the dot-dashed lines and 
B| = 8 for the short-dashed curves 



We have analyzed the impact of the various DAs 
of the T^'-meson on the VF. The quark component of 
the VF is stable for different values of G [0, 8] 
(Fig. [H panel (a)). In contrast, the gluon component 
of the VF demonstrates rapid growth with i3| (Fig. 
|6)d). As a result, due to different signs of the quark 
and gluon components of the VF, the total vertex 
function Fjj' gg. {Q^ , uj — ±1,77) for _B| ^ runs be- 
low the asymptotic one (Fig. 7a). For comparison 
the predictions derived in the standard HSA are also 
shown (Fig. 7b) . The quantitative difference between 
the corresponding curves is clear. 

The dependence of the full VF on the asymmetry 
parameter co is depicted in Fig. [51 In the calculations 
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Fig. 7. The fuU VF at tj = ±1 computed using the RC 
method (a) and the standard HSA (b). The correspon- 
dence between the depicted lines and the parameter Bf 
is the same as in Fig. [S] (the short-dashed lines are not 
shown) 
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Q) 

Fig. 8. The full VF obtained by employing the RC 
method (the solid lines) and the standard HSA (the 
dashed lines) as a function of lo. For the lines 1 the pa- 

1 = 4, and for the 
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rameter is B| = 0; for the lines 2 
lines 3 - Sf 



the framework of the standard HSA. For the sake of 
simplicity, let us consider the ?7'-meson asymptotic 
DA. In this case, the standard pQCD prediction for 
the VF is given by ([50)1 . Now, we want to approximate 
the RC prediction in the form 



F^%.^.(Q^w,77) = - ""^3gr ' K2{u;,rj) (35) 

b c 



X2^i(1.2;4;r2)(a+ — + ^ 



From the general formula ([35]), for the ry'-meson on- 
shell gluon transition we get 



the r?'-meson DAs with various values of Bf are em- 
ployed. As it has been noted above, the gluon compo- 
nent of the VF is identically equal to zero at u — 0, 
and the quark component, as a function of to, demon- 
strates stability against variations of B|. Therefore, 
it is easy to understand the features of the full VF 
as a function of uj. Really, in the region \uj\ < 0.3 the 
difference between the VFs corresponding to differ- 
ent i?f is small; it becomes essential for > 0.8. 
But, owing to the power corrections, mainly to the 
quark component, even for < 0.3 the Borel re- 
summed full VF significantly exceeds the standard 
pQCD result. 

For phenomenological applications it is useful to 
parameterize the VF using some simple expressions. 
We start from the expressions of the VF obtained in 



F", . 
V gg 



±1,77) 



(36) 



1 + 77 



The fitting procedure gives the following values of the 
parameters: 



1.7837, 6- 0.8228, 



-1.022. 



In the case of the gluons with equal virtualities, 
we find 



V g g 



3g2 



Koiv) (37) 
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Fig. 9. The approximations to the full VF. For the upper 
curves ui — 1, for the lower ones ui — 0. The solid lines 
are RC predictions, the dashed ones the corresponding 
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X2Fi(l,2;4;ro) (^a + ^2 + q4 + ■ ' ' 

with 



We have extended the results obtained in \22} for 
the massless ry'-meson virtual gluon transition VF. 
It has been shown that effects generated by the i]'- 
meson mass term considerably change the predictions 
obtained in [22 : they suppress the absolute values 
of the quark and gluon components of the VF, and 
modify their behavior as functions of the asymmetry 
parameter lu. This modification, in the case of the 
gluon component, has not only a quantitative, but 
also a qualitative character: thus, at lu — the gluon 
component of the VF vanishes identically. As a re- 
sult, mass effects change the dependence of the full 
VF on the total gluon virtuality and asymmetry 
parameter lo. 

The numerical analysis presented shows that power 
corrections considerably enhance the standard pQCD 
predictions for the VF in the explored region 
1 GeV^ < Q2 < 25 GeV^, though other sources, may 
also give rise to power corrections. As an important 
consistency check, we have proven that the results ob- 
tained within the RC method in the asymptotic limit 

— > oo reproduce the standard pQCD predictions 
for the vertex function. 



a - 1.5485, 5- 2.3361, c~ -1.627. 

The corresponding results are shown in Fig. (5] As is 
seen, p7l) leads to an almost perfect approximation 
of the original result, whereas the expression ([36|) de- 
scribes the exact prediction, demonstrating, never- 
theless, some deviations. 

It is known, that the RC method produces higher- 
twist ambiguities. For the massless 77'-meson gluon 
transition VF they were estimated in [22] and found 
to lie within ±15% of the original results. Because 
such modifications cannot change our principal con- 
clusions, we do not concentrate on these questions 
here. 



6 Concluding remarks 

In this paper we have evaluated soft end-point (power- 
suppressed) corrections to the space-like ry'-meson vir- 
tual gluon transition VF by including the jy'-meson 
mass effects. To this end, we have employed the stan- 
dard HSA and RC method in conjunction with the 
IR renormalon calculus. In the calculations, both the 
quark and the gluon components of the ?7'-meson DA 
have been taken into account. We have modelled the 
DAs by retaining in the general expressions ([6]) and 
([7|) only the first non-asymptotic terms. 
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